analyte_plot <- function(df_i,analyte_i,regime_i) {
df_i %>%
filter(Analyte_Name==!! analyte_i & Puffing_Regime==!! regime_i) %>%
ggplot(aes(y= Unformatted_Result, x=Test_Article_ID, color=Puff_Collections)) +
geom_boxplot(alpha=1)+
scale_colour_brewer(palette = "Set1")+
labs(
y = "mass analyte/50 puff collection",
x = "Sample Batch")+
ggtitle(paste0("Aerosol ",analyte_i," Level, 0 Month, ", regime_i, " Puffing Regime"))+
theme(plot.title = element_text(size = 12, face = "bold"),
axis.text = element_text(size = 6),
axis.title = element_text(size = 10))+
scale_y_continuous(breaks=number_ticks(10))+
theme(legend.position = "right",
axis.text.x = element_text(angle = 30, hjust=1))+
facet_wrap(~Formulation_Name, scales = "free")
}
#e,g. analyte_plot(MV_v1_PMTA_T0, 'Nicotine','55/3/30')
unique_analytes <- as.vector(unique(MV_v1_PMTA_T0$Analyte_Name))
unique_analytes_plot <- unique_analytes[! unique_analytes %in% c("Impinger Solution Volume", "Number of Impingers")]
unique_regimes <- as.vector(unique(MV_v1_PMTA_T0$Puffing_Regime))
unique_combquery <- expand.grid(unique_analytes_plot, unique_regimes,stringsAsFactors = FALSE)
for (i in seq_len(nrow(unique_combquery))) {
print(analyte_plot(MV_v1_PMTA_T0, paste0(unique_combquery[i,1]), paste0(unique_combquery[i,2])))
}






















































































